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Factorization analysis for the fragmentation functions of 
hadrons containing a heavy quark 



a = Q,Q,q,q,g). By integrating out pair production of virtual and real heavy quarks, 
the fragmentation functions are matched onto a single nonperturbative function describ- 
ing the fragmentation of the heavy quark Q into the hadron H in "partially quenched" 
QCD. All calculable, short-distance dependence on x is extracted in this step. For 
x — ► 1, the remaining fragmentation function can be matched further onto a universal 
function defined in heavy-quark effective theory in order to factor off its residual de- 
pendence on the heavy-quark mass. By solving the evolution equation in the effective 
theory analytically, large logarithms of the ratio n/mq are resummed to all orders in 
perturbation theory. Connections with existing approaches to heavy-quark fragmenta- 
tion are discussed. In particular, it is shown that previous attempts to extract ln n (l — x) 
terms from the fragmentation function Dq/jj(x,171q, fj,) are incompatible with a proper 
separation of short- and long-distance effects. 
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1 Introduction 



The production of heavy quarks (bottom or charm) in high-energy collisions and their sub- 
sequent fragmentation into heavy hadrons is an important process in particle physics. For 
instance, a light Standard Model Higgs boson would preferably decay into a pair of bb quarks, 
which then fragment into B mesons or heavier beauty states. The direct production of heavy 
quarks would be an important background process, so a precise calculation of the heavy-quark 
production and fragmentation processes becomes mandatory. Since the heavy-quark mass 
provides a natural infrared cutoff, the calculation of the cross section can, to a large extent, 
be performed using perturbative QCD. 

Since the center-of-mass energy y/s of modern particle colliders is much larger than the 
bottom- or charm-quark masses, the total inclusive heavy-quark production cross section is 
well approximated by neglecting the masses mq of the heavy quarks. However, less inclusive 
processes such as the differential energy distribution of heavy hadrons produced in e + e~ or 
pp collisions are sensitive to logarithms of the ratio s/rriq. The presence of these logarithms 
threatens the convergence of perturbative expansions. Fortunately, it is possible to show 
that up to power corrections in rriq/s the cross section factorizes into mass-independent hard 
partonic cross sections convoluted with fragmentation functions, which describe the probability 
that a parton fragments into a particular hadron containing the heavy quark. 

Consider, for concreteness, the production of a heavy meson H containing a heavy quark 
Q via the decay of a vector boson V = 7, Z produced in e + e~ annihilation, 

e + e- ^V(q)^H(p H )+X, (1) 

where s = q 2 . Working in the center-of-mass frame, the energy of the heavy hadron can be 
measured using the scaling variable 

2Ph ■ q cms 2E H 
x = — = — r . 2 

We will be interested in a situation where s ^> trig. At leading power in tjiq/s, but to all 
orders in perturbation theory, the differential production cross section can be factorized as [1] 
(see also [2]) 

= E^( I »v / s^)®- D «/fi( I . m O,M), (3) 

where D a /H(x,mQ, n) gives the probability for a parton a to fragment into a heavy meson H 
carrying a fraction x of the parton's momentum, and da a /dx is the cross section for producing 
a massless parton a with energy fraction x defined in analogy with (2). The symbol <g> denotes 
the convolution 

A(x) ® B(x) = ( 1 dz l [ 1 dz 2 A(z 1 )B(z 2 )5(z-z 1 z 2 )= f 1 —A(z)B(-). (4) 

JO JO Jx Z K Z> 

The factorization formula (3) separates the dependence on the heavy-quark mass uiq from 
the dependence on the center-of-mass energy, which is contained in the partonic cross sec- 
tions. These hard cross sections can be calculated in the massless approximation. At lead- 
ing power all mass dependence resides in the fragmentation functions, which are universal 
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(i.e., process-independent) nonperturbative quantities. Like the parton distribution functions, 
fragmentation functions must be extracted from one set of processes and can then be used to 
calculate the rates for other processes. Large logarithms of the ratio s/itiq can be resummed 
to all orders in perturbation theory by evolving the fragmentation functions from a low scale 
/i ~ rriQ up to a high scale /i ~ \/s by solving the DGLAP evolution equations 



d 



dlnfi'- 



D a/H (x, m Q , fi) = Pa^b(x, fi) <g> D b/H (x, rriQ, ji) , 



(5) 



where P a ^b are the time-like Altarelli-Parisi splitting functions [3, 4, 5, 6]. Alternatively, the 
partonic cross section can be evolved from the high scale down to a scale of order toq. At 
leading order in perturbation theory (but not beyond), the functions P a ^b are related to the 
space-like splitting functions Pb+- a governing the evolution of the parton distribution functions 
by the Gribov-Lipatov reciprocity relation [7]. 

In contrast to the case of parton distribution functions, the moments of the fragmentation 
functions cannot be related to local operator matrix elements. Even the normalization of 
the leading fragmentation function Dq/h is unknown, since a heavy quark Q can fragment 
into different heavy hadron states H. In the literature on heavy-quark fragmentation one 
frequently encounters the notion of "perturbative fragmentation functions" D a /q [1], which 
describe the probabilities for partons a to fragment into an on-shell heavy quark Q. These 
functions are finite in perturbative QCD, since the heavy-quark mass, or more precisely the 
scale mq{l — x), provides an infrared cutoff. They are relevant for the discussion of inclusive 
heavy-quark production, where one sums over all possible hadron states H containing the 
heavy quark Q. Indeed, quark-hadron duality suggests that 



D a/H (x,m Q ,ti) = D a/Q (x,m Q ,fi) + 0[ 

H 



A 



QCD 



\m Q (l -x) 



(6) 



where the sum is over all hadron states containing the heavy quark Q. However, such a relation 
can be expected to hold only if x is not too close to 1, such that the scale mg(l — x) 3> Aq C d is 
in the short-distance regime. The one-loop results for the perturbative fragmentation functions 
were obtained a long time ago [1]. They read 



D Q/Q (x, m Q , fj) =6(l-x) + 



C F a s 
2vr 



1 + x 2 
l-x 



In 



m 2 Q {l-xy 



- 1 



+ 0(a 2 s ) , 



D g/Q (x,m Q ,fi) = 



T F a s 



x 



+ {l-xf 



ln-^ + O^ 2 ), 

m Q 



(7) 



while D a /Q(x,mQ, fi) = 0{a 2 ) for a ^ Q,g. Throughout our paper (as is the case in most of 
the literature on heavy-quark fragmentation) itlq denotes the pole mass of the heavy quark. 
If desired, this parameter can be expressed in terms of a short-distance mass parameter using 
perturbation theory. The two-loop expressions for D a /Q(x,rriQ, /i) have been calculated in 
[8,9]. 

The strong dynamics binding the heavy quark inside a heavy meson implies that there 
exist important nonperturbative effects leading to a prominent peak of the fragmentation 
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function Dq/ h near the endpoint, 1 — x pcak = 0(Aqcd/"^q)- For example, the popular model 
of Peterson et al. [10] 

, . Nx(l -x) 2 
Dq/h{x) = \ } (8) 

[(1 - X) 2 + € Q X\ 

exhibits a peak at rc pca k ~ 1 — y/eQ, where €q ~ (Aqcd/^q) 2 - It is therefore interesting 
to study the heavy-quark fragmentation process (1) in the kinematic region x — > 1, where 
(1 — x) is treated as a small parameter. In this region the partonic cross sections da a /dx 
contain large logarithms of (1 — x), which need to be resummed to all orders in perturbation 
theory. The traditional approach to resumming these logarithms proceeds via Mellin moment 
space [11, 12]. It has first been applied to the case of the fragmentation process in [13]. 2 
More transparently, it is possible to formulate a factorization theorem for the partonic cross 
sections using methods of effective field theory and to perform the resummation of large 
logarithms directly in x space, applying a novel approach developed recently in [14, 15]. This 
will be discussed in a forthcoming paper [16]. In the x — > 1 limit, the fragmentation functions 
are sensitive to contributions associated with three different mass scales: ttlq, mg(l — x), 
and Aqcd- The dominant contributions to the cross section come from the region where 
7Bq(1 — x) ~ Aqcd, so that it is natural to assign the power counting (1 — x) ~ Aqcd/^q- 
The goal of the present work is to systematically factorize short- and long-distance effects for 
the fragmentation functions D a /H(x,rriQ, fi). The leading singular terms for x — > 1 reside in 
the fragmentation function Dq/h{x, mg, fi) describing the conversion of a heavy quark Q into 
a heavy hadron H containing that quark. 

The problem of how to systematically incorporate and describe nonperturbative bound- 
state effects into the heavy-quark fragmentation function has been addressed previously by 
several authors. The most popular approach is to factorize the fragmentation function into 
perturbative and nonperturbative components: Dq/h — ^q/h ® -^q/h [17, 18, 19] (for a 
recent analysis using this scheme, see e.g. [20]). The first component is identified with the 
perturbative fragmentation function, Dq'J^x, tuq, /i) = Dq/q(x, tuq, /i), while for the nonper- 
turbative component a model such as (8) is adopted. An important observation we make in 
the present paper is that such an ansatz is incompatible with a proper factorization of short- 
and long-distance contributions. For x — > 1 the perturbative fragmentation function itself con- 
tains long-distance contributions from momentum scales of order wq(1 — x) ~ Aqcd, which is 
evident from the form of the logarithmic term in the first relation in (7). Such logarithms are 
not controllable in perturbation theory and hence must be factorized into the nonperturbative 
function D^ H . Previous attempts to resum the ln n (l — x) terms in the fragmentation func- 
tions have thus led to incorrect formulae for the perturbative Sudakov factors in the x — > 1 
region. To derive the correct form of the factorization formula for the fragmentation functions 
is the main objective of our work. 

As an aside, we mention that the situation encountered here resembles that in the theo- 
retical analysis of inclusive I?-meson decays into light hadrons in the "shape function region" . 
There as well the systematic factorization of short- and long-distance contributions using ef- 
fective field-theory methods [21, 22] has helped to identify some long-distance contributions, 
which had previously been included in the perturbative decay rate [23] . The proper subtraction 



2 See also [1] for a discussion of the complications arising in the x — * 1 region. 
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of these terms had a significant impact on phenomenology [24]. 

Alternative schemes attempt to model the hadronization process by introducing nonper- 
turbative effects into the (resummed) perturbative fragmentation functions themselves, for 
instance by using a renormalon analysis [25, 26] or by introducing an effective, infrared-safe 
low-energy QCD coupling constant [27, 28]. A different approach was introduced in [29], 
where the long-distance effects in the fragmentation process were described using heavy-quark 
effective theory (HQET). A subsequent work followed up on this idea and aimed at model- 
ing the bound-state effects in the heavy hadron in terms of hadronic form factors [30]. The 
compatibility of the HQET-inspired approaches and those based on perturbative QCD was 
studied in [31]. Only much later the first systematic matching between the two approaches 
was performed, and the evolution equation for the HQET fragmentation function was derived 
[32]. In the same paper, a simple relation with the evolution equation for the 5-meson shape 
function [33, 34] was discovered. In fact, there is a close analogy between the factorization 
analyses for the heavy-quark fragmentation process and that of inclusive 5-meson decays into 
light particles, such as B — > X s ^ and B — > X u lu [21, 22, 24]. Using the methods developed in 
these papers, it would be possible to describe power corrections to the fragmentation functions 
D a m in the x — > 1 limit in terms of subleading fragmentation functions in HQET, which are 
analogous to the subleading 5-meson shape functions studied in [35, 36, 37]. Such an analysis 
is, however, beyond the scope of the present work, and we argue that it would be of little 
phenomenological relevance. 

The remainder of this paper is organized as follows: In the following section we perform a 
first factorization step at the hard scale uiq by integrating out real and virtual heavy-quark 
pairs. This results in a "partially quenched" version of QCD, in which the heavy quark that 
fragments into the heavy hadron is still present, but no pair production of heavy quarks is 
allowed. In this first matching step the different fragmentation functions D a /u are matched 
onto a single, nonperturbative function describing the fragmentation of the valence heavy 
quark Q into the heavy hadron H . This matching step can be performed irrespective of the 
value of x. From a practical point of view, it is the most important result of our analysis, 
since at this stage all calculable dependence on x has been extracted from the fragmentation 
functions. In Section 3 we then consider the limit x — > 1 and perform a second factorization 
step at the hard scale m<g, in which "partially quenched" QCD is matched onto HQET. This 
is appropriate whenever mg(l — x) <C uiq, and it is a natural step if tuq^I — x) ~ Aqqd, as is 
the case for the region near the peak of the heavy-quark fragmentation function. We rederive 
the relation between the perturbative fragmentation function and the perturbative 5-meson 
shape function noted in [32] and use it to extract the two-loop matching coefficient relating 
the fragmentation functions defined in HQET and "partially quenched" QCD. We also present 
the exact analytic solution of the evolution equation obeyed by the fragmentation function, 
which is analogous to the solution of the evolution equation for the shape function obtained in 
[38] . In Section 4 we illustrate the phenomenological implications of our results with the help 
of a simple model for the fragmentation function in HQET. We study the scale dependence 
of the fragmentation function Dq/h and derive a heavy-quark symmetry relation between the 
fragmentation functions of B and D mesons. Section 5 contains a summary and conclusions. 
Some technical details of our analysis are relegated to three appendices. 
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Figure 1: Examples of fragmentation processes converting a parton a into a heavy hadron H. 
Heavy quarks are denoted by thick lines, while the heavy hadron is drawn as a double line. 
Only the first process is allowed in "partially quenched" QCD. 

2 Decoupling of real and virtual heavy-quark pairs 

The fragmentation functions D a /H(x,rriQ, fi) receive contributions associated with different 
mass scales, such as thq, mg(l — x), and Aqcd- Our goal is to factorize short-distance effects 
related to the hard scale tuq from contributions related to lower scales. It turns out to be 
useful to integrate out the heavy-quark scale tuq in two steps. In the first step, QCD with 
rif = rii + 1 active flavors, in which the factorization formula (3) has been derived, is matched 
onto a theory with ni active flavors, which no longer contains QQ pair production or heavy- 
quark loops. As always, light flavors are treated as massless. In the MS scheme, the matching 
relation between the running coupling constants a s (fj,) in the two theories reads 



(rii+l) 



(m) 



a 



(ni) 



2 m a\ 
1 H — Tp — 
3 



2tt 



+ 



(9) 



We would like to write an analogous matching relation for the fragmentation functions D a m. 
In the theory where the heavy quark Q is among the active flavors, the fragmentation process 
can be initiated by any parton a = Q,Q,q,q, g, as illustrated in Figure 1. However, in all cases 
but a = Q this will involve the production of a QQ pair, which is forbidden in the n/-flavor 
theory. The most general matching relation for the fragmentation functions defined in the two 
theories therefore reads 



DT£ l \x,m Q ,n) 



'a/H v*>'"«>fv - C a/Q(x,m Q ,ii)®D i Q} ) H (x,mQ,iJ l ). (10) 
The meaning of this relation is obvious: the functions C a /Q describe the fragmentation of a 
massless parton a into a heavy quark Q, while Dq) h describes the fragmentation of this heavy 
quark into the heavy hadron H. Only the latter process is sensitive to long-distance dynamics. 
The functions C a /Q can be extracted from the perturbative expressions for the fragmenta- 



tion functions derived in [8, 9]. We obtain 
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Cq/q(x, m Q , n) = 6(l-x) + cq(x) ■ 
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a = Q,q,q. 



(11) 
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To the extent they are known, the expressions for the two-loop coefficients c a (x) are collected in 
Appendix A. These coefficients receive contributions from two sources: terms in the fragmen- 
tation functions of the (rii + l)-flavor theory that involve a heavy-quark loop or heavy-quark 
pair production, and terms resulting from the replacement of the effective coupling constant 
in the one-loop fragmentation functions according to (9). The leading terms in the x — > 1 
limit reside in cq and are given by 



1 2 /i 2 (I 2tt 2 \ /i 2 3139 7T 2 2 ' 



c Q (x) = C F T F - x) 

1 /2, 2 pi 2 20 , n z 56'. . . ,. 

+ -— -ln^-- T ln^ + - K subleadmg terms. (12) 



The remaining two-loop matching corrections in (11) are regular for x — > 1. It follows from 
(5) that the matching coefficients C a /Q and the remaining fragmentation function Dq 1 ^ h obey 
the evolution equations 



An rUl^ 

gWg _ \- p( n < +1 ) r r (9> p( n; ) Q/g _ p(n.) ^ n (n,) /,o\ 

^2-L^6 ®C 6 /g-C7 a/0 ®P^ , "cHn/i 2- ~~ Q^Q Q/h ■ (13) 

The first relation can be checked using the explicit results for the Wilson coefficients C a /g 
given above along with the known expressions for the Altarelli-Parisi splitting functions. 

The short-distance coefficients C a /Q can be absorbed into the definition of a partonic 
heavy-quark production cross section defined in the nrflavor theory: 

da { n l) da^+V 

—^-{x,y/s,m Q ,n) = ^ — {x,y/s,n)®C a/Q (x,m Q ,n) . (14) 

The differential heavy-hadron production cross section can then be written in the simple form 

da da^ 

~~dx~ = "cfo - ^' mQ ' ^ ® ^S^' m Q' V) ■ ( 15 ) 

The advantage is that now only a single fragmentation function remains, which contains all 
nonperturbative information. The cross section da^ jdx contains large logarithms, which can 
be resummed by evolving the partonic cross sections da^ 11 ^ /dx from a high scale down to a 
scale n ~ rriQ using the DGLAP evolution equations. 

From a practical point of view, eqs. (10), (14), and (15) are the most important re- 
sults of our work. While in the following section we will factorize the fragmentation func- 
tion DQ l ^ H (x,mQ, n) further by factorizing off its remaining dependence on the heavy-quark 
mass, all calculable short- distance dependence on x has already been extracted in the stage 
above. For phenomenological applications it would be appropriate to treat the function 
Dq^ h (x,i71q, h) as a hadronic quantity that is extracted from a fit to experimental data. The 
discussion of the next section becomes relevant only when one tries to relate fragmentation 
functions describing different processes (see Section 4). 
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At this point a comment is in order about a technical detail in the matching calculation 
which, at first sight, appears problematic. When integrating out heavy-quark pairs we com- 
pute short-distance effects associated with the hard scale vtiq. However, the resulting Wilson 
coefficient Cq in (12) contains a term proportional to 1/(1 — x)+, which for x — > 1 generates a 
large logarithm. Indeed, an integral of the original fragmentation function over the endpoint 
region with width (1 — x ) ~ Aq CD /mQ appears to generate an infrared-sensitive logarithm, 



which is not accounted for by the fragmentation function in the low-energy effective theory. 
But this would be in conflict with the general fact that a low-energy effective theory must 
be identical to the full theory in the infrared domain. The resolution of the puzzle is that 
the ln(l — Xo) term in the relation above should not be interpreted as an infrared-sensitive 
logarithm arising from the ratio of the two scales rriQ(l — x) ~ Aqcd and itlq, but as an 
ultraviolet logarithm arising from the presence of the two large scales s and s(l — x) in the 
fragmentation process (1), where s(l — x) is the invariant mass squared of the final state X, 
which we assume to be much larger than tuq. In fact, these are the only large logarithms that 
can be extracted from the fragmentation functions in a meaningful way. In a forthcoming paper 
we will show that the /z-dependent terms in (12) combine with other terms in the resummed 
expression for the partonic cross section dag/dx in such a way that the running coupling 
a s (fx) is converted from the (rii + l)-flavor to the n r flavor theory [16]. The constant pieces 
give rise to matching corrections, which are analogous to the matching correction creating a 
discontinuity in the running coupling constant at three-loop order. 

From now on, we focus on the fragmentation function -Dq/# defined in the theory with n\ 
active flavors. In the discussion below, quantities without superscript will always be defined 
in the theory with n\ massless flavors, i.e., we use a s = and Dq/h = Dq L ^ h . 

3 Factorization for x — ► 1 

The step of integrating out virtual and real heavy-quark pairs discussed in the previous section 
leads to the factorization formula (10), in which the different fragmentation functions of full 
QCD are related to a single nonperturbative function Dq/ H defined in the "partially quenched" 
n;-flavor theory. This first matching step can be performed irrespective of the value of x. In 
particular, relation (15) for the fragmentation cross section is exact. 

We now return to the specific case of the endpoint region, where (1 — x) ~ Aqcd/^q is 
treated as a small parameter. We wish to factorize effects associated with the short- distance 
scale rriQ, which are still present in the n r flavor theory, from long-distance effects associated 
with the scales m<g(l — x) and Aq CD . HQET is the appropriate effective theory for describing 
these nonperturbative effects [39]. Matching "partially quenched" n r flavor QCD onto HQET 



['dxD^&mQ,^ 3 C F T F 




) 



(16) 
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we obtain a matching relation of the form 

D Q/H {x,m Q ,fi) -> C D {m Q ,fi) \ S Q/H {u, fi) + — ]T Cf } (m Q , /i) sf } (^, ft) + ...I , (17) 

where u> = O(Aqcd) is a momentum variable in the effective theory yet to be defined. The 
Wilson coefficients Co and depend on the heavy-quark mass logarithmically through the 
running coupling a s (mQ). The functions Sq/h and S\ are effective, leading and subleading 
fragmentation functions defined in HQET (see below). The contributions of the subleading 
functions are suppressed with respect to the leading term by powers of Aqcd/^q or (1 — x). 
For most of this section we will be concerned with the leading term. 

We employ the conventional definition of the (bare) fragmentation function (defined in 
d = 4 — 2e space-time dimensions) of a heavy quark Q transforming into a heavy hadron H 
with light-cone momentum fraction x, given by [5] 



x' 



d-3 



D Q/H (x,m Q ,ri = ^J ^e^^^(0|(^ Q ^)^M|//Mx>(//Wx|(C/„ r ^)L(0)|0). 

(18) 

Here n M is a light-like vector, a, f3 are Dirac indices, and i is a color index. Longitudinal boost 
invariance implies that the result should be unchanged under rescalings of the light-like vector 
n. It follows that the right-hand side of (18) can only depend on the ratio x = n-pn / n-p G [0,1]. 
The expression above is rendered gauge invariant by the insertion of Wilson lines. We define 

U n (y) = Pexp (ig J ds n ■ A(y + sn)^j , U ] n {y) = P exp (-ig J ds n ■ A(y + sn)^j , (19) 

where the symbol P means path-ordering, such that gluon fields are ordered from left to right 
in the order of decreasing s values, and P means the opposite ordering. In (18) we need the 
related objects U*(y) and U^(y), which are obtained from U^{y) and U n (y) by transposition. 
This reverses the ordering and replaces the gluon fields by the transposed fields A T = A a t^, 
where t a are the hermitian generators of color SU (N c ) . The transposed Wilson lines appear 
because in (18) the Wilson lines are located on the "wrong" sides of the quark fields. 

The fragmentation function receives contributions from both hard and soft interactions. 
In order to separate short- and long-distance effects, we match Dq/ h onto a corresponding 
function defined in HQET. Following [29], we define 

S q/h {u,vl) = ^J rfte*"'E(0|(^^)LM|^oo(^)x s )(^oo(^)x s |(^)L(0)|0) , (20) 

where v is the four- velocity of the heavy hadron, h v (y) is the two-component heavy-quark 
field in HQET [39], and we have used that between these fields iji can be replaced with 
n ■ v. The hadron states in HQET are taken to have a mass-independent normaliza- 
tion (HoolhyhvlHoo) = 1 instead of the conventional normalization to 2Mh employed in (18). 
In the effective theory we sum over soft states Xs (i-e., states having momenta much less than 
mg), and the QCD Wilson lines U n , are replaced by soft Wilson lines S n , defined in 
analogy with (19), but with QCD gluon fields replaced with soft gluon fields. Note that k 
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corresponds to the residual momentum of the heavy quark, defined as p = tuqv + k. Boost 
invariance implies that the right-hand side must be a function of uj = n ■ k/n ■ v G [A, oo[, 
where A = Mu — itlq is the residual mass of the heavy hadron state in HQET, and the condi- 
tion uj > A follows since the heavy quark must have sufficient energy to produce the heavier 
hadron state H. The above definitions imply the relation uj = Mh/x — rriQ, and indeed one 
of the main points of [29] was to argue that the fragmentation function should be considered 
as a function of this particular combination of x, tuq and Mh (see also [31]). We find it more 
convenient to write the fragmentation function instead as a function of the variable uj = uj — A, 
which takes values between and oo and obeys the simpler relations 

a = *=^. ( 21) 

x M H + uj 

After the decoupling of heavy-quark pairs discussed in the previous section, the matching 
relation between the fragmentation functions in QCD and HQET is local in x and u>, as shown 
in (17). Using the usual definition of the HQET spinor h v [39] one readily obtains the tree- 
level relation Dq/ h (x, ttiq, //) = Mh x Sq/h(&, /-*) + •••, where the factor M# results from the 
different normalization of the hadron states in (18) and (20), and the factor x is due to the 
x d-3 p re f ac tor in (18). Beyond tree level, we write the matching relation in the differential 
form 

dx 

Dq/ h (x, tuq, fj) — = g(x) C D {m Q , //) S q/h (uj, fi) du + power corrections. (22) 

Ob 

Any deviation of the function g(x) from 1 is formally a power correction, so instead of the form 
g(x) = x 2 suggested by the tree-level matching relation given above it would be consistent to 
set g(x) = 1 or use any other smooth function obeying the constraint g(l) = 1. We will see in 
Section 3.4 that a particularly convenient choice is g(x) = (1 + x 2 )/2. 

General arguments analogous to those presented in [22] suggest that the HQET fragmen- 
tation fragmentation vanishes for uj — > 0, which implies that the QCD fragmentation function 
vanishes for x — > 1. Note, however, that in general there is no reason why the fragmentation 
function should vanish at x — 0, even though this is built into most phenomenological param- 
eterizations. Indeed, the perturbative results in (7) and the corresponding two-loop results 
in [8, 9] suggest that the fragmentation functions tend to constants modulo logarithms for 
small x. 



3.1 Matching at one-loop order 

The hard matching coefficient Cd can be calculated in perturbation theory by computing the 
fragmentation functions in the two theories and requiring that the two sides in (22) be identical 
up to power-suppressed terms. Since the effective theory (HQET) exactly matches the full 
theory ("partially quenched" QCD with n\ light flavors) in the infrared, the matching can be 
performed using on-shell quark and gluon states. Examples of one-loop diagrams contributing 
to the perturbative fragmentation functions in the two theories are shown in Figure 2. 

Replacing the heavy hadron by an on-shell heavy quark, we find that at one-loop order in 
n;-flavor QCD the renormalized fragmentation function Dq/q in the MS scheme is given by the 
first formula in (7). At the parton level we do not distinguish between the heavy-quark mass 
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Figure 2: Examples of one-loop diagrams contributing to the perturbative fragmentation 
functions. The grey dots represent the Wilson lines in expressions (18) and (20). 

and the mass of the hadron containing the heavy quark (i.e., we set M# = tjiq). Performing 
the corresponding on-shell calculation in HQET, we obtain 



Sq/q(&,(J>) =S(u)[l 



CpCi s 7T \ CfO-s 



I (i + 21n* 



(23) 



7T 24 / 7T 

where the star distributions are defined as [22, 23] 

fa m (twm) M = fa m-m h . £ + /a ln „ +1 5 . (24) 

jo y w / ^ jo u; /in+1 n 

This generalization of the usual plus distributions is required, because in the effective theory 
there is no reasonable way of defining a dimensionless variable x G [0, 1]. 

The expressions (7) and (23) must be understood in the sense of distributions. In order to 
extract the short- distance matching coefficient Co, we integrate these perturbative fragmen- 
tation functions with an arbitrary test function over an interval x G [mg/(mg + ujq), 1] and 
Cj G [0, wo], respectively, where ojq <C tuq. We then require that the results agree to leading 
power in uq/uiq. This leads to 

C D(mQ , ,) = ! + ^ (I ,„ 2 £l + I ,„ £l + 2 + g + (qS . (25) 

This result agrees with a corresponding expression derived in [32]. We stress again that the 
matching coefficient is free of m(l— x) terms. All those terms in the perturbative fragmentation 
function in (7) are of long-distance origin and are still contained in the HQET fragmentation 
function S Q / H (u,fj)- 



3.2 Higher-order matching 

It is an intriguing observation that the expression (23) for the perturbative HQET fragmen- 
tation function can be related to the corresponding perturbative expression for the 5-meson 
shape function. The shape function is the HQET analog of the heavy-quark parton distribu- 
tion function of the hadron H [33, 34, 40]. In analogy with (20), it is defined as 

77 * 1) C — 

S H (u>,n) = — ydte ifc ^(^ oo (7;)|(^S' n )(0)(5't/ lt; )(in)|^ oo (7;)>. (26) 
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In this case uj = n ■ k/n ■ v is restricted to take values between — oo and A, since the residual 
light-cone energy of the heavy quark can at most equal the residual mass A of the heavy 
hadron. 3 We thus define Co = A — uj e [0,oo[. When this is done, our one-loop expression 
in (23) coincides with the one-loop expression for the perturbative shape function given in 
eq. (33) of [22]. 

The agreement between the perturbative fragmentation and shape functions in HQET, 
when expressed in terms of the proper variables, is not restricted to one-loop order. Rather, 
it is an all-order identity between the two functions, which follows from their relation to 
vacuum expectation values of Wilson loops. This has also been observed in [32]. In order 
to derive this relation we make use of field redefinitions, which decouple the interactions of 
soft gluons from the effective heavy-quark fields. For the case of the shape function, we use 
h v (y) = S v (y) h^(y) (see, e.g., [41]), where the soft Wilson line S v is defined in analogy with 
S n , but with the light-like vector n replaced with the time-like vector v. For the case of the 
fragmentation function we need the Wilson lines on the "wrong" side of the heavy-quark field, 
so we use instead h v (y) = h^\y) £F v (y). The new fields are "sterile" and do not couple 
to any other fields in the theory. The corresponding field operators simply create or annihilate 
the heavy quarks in the external hadron states, giving rise to free two-component Dirac spinors 
u v satisfying tfu v = u v and u v u v = 1. We define 

h^lHM) = \<f H (Av)) , (27) 



where the state \q l H (Av)) represents the light degrees of freedom inside the hadron H with color 
index %. This state carries momentum Av, which is what remains from the hadron momentum 
Mhv when the heavy quark is removed. We may therefore write 

S h (u,,vl) = ^Jdte* K ^>*ti H {kv)\ [(SlS n )(0)(SlS v )(tnf \^ H (Av)) . (28) 

After the decoupling transformation, the corresponding expression for the fragmentation func- 
tion reads 

Sq,h@,h) = ^nJ <fte i(A ^ W £<0|(^S^ 

C Xs 

(29) 

In a partonic picture the spectator quark does not participate in the short-distance process, 
and to obtain the perturbative shape and fragmentation functions it can simply be dropped 
from the above matrix elements. At the same time the parameter A is set to zero, thereby 
removing the contribution of the spectator quark to the overall momentum balance. This 
yields for the perturbative shape function 

S Q (^,/i) = ^/^e-™'^(0|Tr(5X)(0)(^)(tn)|0). (30) 



3 Note that in the original variable w — n-k/n-v the range of support of the shape function is the complement 
of the range of support for the HQET fragmentation function. 
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In the case of the perturbative fragmentation function we can now perform the sum over the 
soft states Xs, since after the removal of the spectator quark this is a complete sum over states 
in the effective theory. The result is 

Sq/q^h) = ^J dte^- nt ^(0\Tr(St v S n )%tn)(SiS^y(0)\0) 
= ±-J ^e^*i-(0|Tr(^ n )M(^-,)(0)|0)* 

= 2^/ ^ e ^' n *^(°l^(^n)(0)(^-,)M|0) 

= hi * e ^^< o i^(^ 5 »)(°)( 5 i 5 «)(-*")i o >= 5 ^^)- ( 31 ) 

The last step follows since the position-space Wilson loop (0| Tr (S_ v S n )(0)(S^S- v )(tn)\0) is 
a function of the variable itn • v, so replacing the sign of v is equivalent to changing the sign 
of t. It follows that the perturbative fragmentation function coincides with the perturbative 
shape function. 

Armed with this result, the short-distance coefficient Co in (22) can be computed to 
two-loop order using existing calculations. The two-loop expression for the fragmentation 
function in the "partially quenched" n/-flavor theory can be extracted from [8] by retaining 
the leading terms in the x — > 1 limit in the expression for the perturbative fragmentation 
function Dq/q(x, tuq, //) and accounting for the first matching step (10) using relations (11) 
and (12). The integral over the two-loop perturbative fragmentation function in HQET is 
obtained from the corresponding integral over the two-loop perturbative shape function given 
in eq. (38) of [42]. Matching the two expressions, we find 

C D (m Q ^) = l + ^-^ - + - + 2 + - 



2 2 12 



+ Cf {^~) [ C F H F + C A H A + T F n f H f ]+0(a s s ), 



(32) 



where L = ln(/i 2 / itlq) , and the expansion coefficients at two-loop order are 

'9 f \ ^ (11 
24 




H f = --- — -\T7 + — ) L -^r--^r--C 3 - (33) 
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If desired, the two matching steps discussed above and in Section 2 can be combined into 
a single matching relation between the original fragmentation functions D a m and the HQET 
fragmentation function Sq/h- In the resulting relation, which is of convolution type, the 
Wilson coefficients C a /Q in (10) are simply multiplied by the coefficient Co in (32). 



3.3 Nonperturbative effects 

While the perturbative fragmentation function in HQET coincides with the perturbative shape 
function to all orders in perturbation theory, such a simple relation will no longer be true 
nonperturbatively. This can be seen from the different ways in which the spectator quark 
appears in expressions (28) and (29). It is also apparent from the fact that the residual 
momentum k of the heavy quark obeys different kinematical restrictions in the two cases. In 
the case of the shape function, n-k varies in the range ] — oo, A] (setting n-v — 1 for simplicity), 
and the equation of motion iv-D h v — of HQET enforces the condition (n-k) = [33, 34, 40]. 
As a result, the first moment of the shape function is given by (u))s H = A = Mh — mg, and 
indeed this relation can be used to define a short- distance, running heavy-quark mass called 
the "shape-function mass" to all orders in perturbation theory [22, 43]. 

For the fragmentation function, on the other hand, the variable n ■ k is always positive 
and varies in the range [A, oo[, so that the average value (n ■ k) > A is a nontrivial hadronic 
parameter. It follows that the first moment of the fragmentation function, (Cj) Sq/h = e H > 
0, is not related to any known parameter and must be regarded as an unknown hadronic 
quantity. This, along with the fact that even the normalization of the fragmentation functions 
is unknown, implies that there is much less theoretical handle on the fragmentation function 
than there is on the shape function. 

Nevertheless, the relation between the perturbative shape and fragmentation functions 
might be employed to motivate some simple models for the fragmentation function. In the 
case of the shape function, a renormalon-inspired model based on an analysis of the divergence 
structure of large-order perturbation theory appears to provide a rather good description of 
different decay spectra in the inclusive processes B — > X s ^ and B — > X u £ v [44, 45]. This model 
starts from the perturbative shape function and introduces nonperturbative effects through a 
single model parameter set by the first moment. This suggests that it might be reasonable 
to use the shape function extracted from the B — > X s 7 photon spectrum as a model for the 
HQET fragmentation function. However, we must account for the fact that the normalizations 
and the first moments of the two functions are different. As a simple model, we thus propose 
the form 

S%JP\u, = N H ST d °\^ H) | a _ h • (34) 

A rather flexible, two-parameter model for the shape function, which is known to provide good 
fits to the experimental data, consists of an exponential times a power of Cj [46]. This leads 
to the form (with b > 0) 
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This function obeys the moment relations 

{' OO P OG POO £^ 

\ F(u) = 1 , du oo F(w) = e H , / do> (c^ 2 - (uj) 2 ) F(u>) = . (36) 
jo j o Jo b 

The parameters e# and b can be tuned to fit the data for the first and second moment of the 
fragmentation function. 

3.4 Residual power corrections 

So far, we have discussed the leading term in the matching relation (22) for the fragmentation 
function, while in general there is an infinite series of power-suppressed contributions as shown 
in (17). Like the leading term, these power corrections can be related to nonperturbative, 
subleading fragmentation functions in HQET. Once a basis of such functions has been 
constructed, their Wilson coefficients can be determined perturbatively by matching 
perturbative expressions for the HQET fragmentation functions with the subleading terms in 
the perturbative fragmentation function Dq/q (as well as more complicated functions such as 
Dq/ 9 q) in nj-flavor "partially quenched" QCD, in analogy with our treatment in Section 3.1. 

Using the relations between plus and star distributions collected in Appendix B, we find 
that at one-loop order in perturbation theory the power-suppressed terms not accounted for 
by the leading contribution in the matching relation (22) are given by 

D Q/Q (x, rriQ, n) - g(x) C D (m Q , /i) S Q/Q (u, fi) x 
C F a s f 2 f l+x 2 \f /i 2 \ Mx) Ins | 

The argument of the logarithm once again indicates that this difference is associated with 
physics at the low scale mq(l — x) ~ Aqcd- A particularly simple expression is obtained if we 
use the special form g(x) = (1 + x 2 )/2 for the weight function on the right-hand side of the 
matching relation (22), since in this case the entire first term vanishes. This particular choice 
will be adopted in our numerical analysis in Section 4. 

To perform the matching onto subleading fragmentation functions is beyond the scope of 
the present work. For practical applications of our results, the explicit inclusion of higher- 
order power corrections is not of much importance due to the fact that even the leading-order 
function Sq/h(<^,p) is not constrained by any useful conditions. Unlike the case of the shape 
function, the moments of the fragmentation function are not related to local HQET operators. 
For practical purposes we may therefore define our model for the function Sq/h{oj, /•*) to include 
all power-suppressed terms in (17) - in any case this model must be tuned to fit experimental 
data. This definition introduces some subleading dependence on the heavy-quark mass mq 
into the fragmentation function, while the original Sq/h(£>,h) was strictly a universal, ijiq- 
independent function. 
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3.5 Renormalizat ion-group evolution 

As discussed earlier, the fragmentation function Dq/ h in the nz-navor QCD theory obeys the 
second evolution equation in (13). The relevant splitting function takes the form [15] 

Pq^q{x, fi) = ^ cusp Y + 7^) 6(1 -x) + subleading terms, (38) 
I 1 _ x ) + 

where the regular terms are subleading for x — > 1 . Here r cusp is the cusp anomalous dimension 
of Wilson loops with light-light segments [47], which is known to three- loop order [48]. The 
quantity 7^ can be identified with a combination of anomalous dimensions of operators defined 
in soft-collinear effective theory and is known to the same order [14, 15, 48]. 

It follows from the discussion of Section 3.2 that the fragmentation function in HQET 
obeys the same integro-differential renormalization-group equation as the shape function. It 
reads 

dlnj2 2 Sq/h{u, h) = - j o du' F(oj, Qj\ \i) S q/h (oj', fi) , (39) 

where to all orders in perturbation theory the kernel has the form [22, 49] 

/ 1 \ M 

F(u, u\ n) = -r cusp (a s ) + 7 s (a s ) 6(u> - u') . (40) 

The anomalous dimension / ~f s is known to two-loop order [32, 42]. Finally, from (38) and (40) 
it follows that the hard matching coefficient obeys the evolution equation 



din// 



^C D (m Q ,fj) 



r cusp (a s ) In -H- + ^(a s ) + 7 s (a s 

m Q 



C D {m Q ,n). (41) 



The renormalization-group equations (39) and (41) can be solved analytically and in closed 
form [38]. The solution for the Wilson coefficient reads 

(m \ 2ar ^ h '^ 

C D (m Q ,fjL) = exp \-2S(fjL h ,fjL) - 2a^(fjL h , fi) ~ 2a 7 s(^ h ,^)l — C D {m Q ,fi h ) , (42) 

\Vh ) 

where ~ toq is a hard matching scale, at which the initial condition ^(mg,///,) is free of 
large logarithms. The evolution functions are given by 4 

S(v,A = - f da-^> f a r (w) = - f da l -^>, (43) 

J (3(a) J P{a') J [3(a) 

a a [y) a s (u) a s (u) 

and similarly for the functions and a 7 s. Here f3 = da s /dla/j, is the QCD /3-function. 
The exponential in (42) accomplishes the resummation of large logarithms to all orders in 
perturbation theory. The perturbative expansions of the anomalous dimensions are collected 
in Appendix C. 



4 The Sudakov exponent S should not be confused with the fragmentation function. 
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The solution of the integro-differential evolution equation (39) has been derived in [38] 
using a technique developed in [50] (see also [22, 51]). One obtains 

S Q/H {u,n) = exp [25(^0,^) + 2a 7 s (^o,^)l / du/ ^S^—^x^z , (44) 

L J 1 {rj) Jo fi '{u — u') 1 n 

where r] = 2ar(/x, fJ-o)- Here /xo serves as a low reference scale still in the perturbative domain, 
at which a model for the HQET fragmentation function is provided. With the help of the 
above relation, the fragmentation function can then be evolved to a higher scale //. It should be 
stressed that the solutions (42) and (44) are formally independent of the two matching scales 
fj, h and fj, . In practice, a scale dependence remains when one truncates the perturbative 
expansions of the evolution functions. 

The renormalizat ion-group improved expressions for the Wilson coefficient and the HQET 
fragmentation function can be combined and simplified to obtain the final result for the frag- 
mentation function Dq/ h in "partially quenched" QCD. We find 



M H g(x) ( /jl y arWh> 

D Q / H (x,mQ,n) = — exp \2S(/j, , fih) + 2a 7 s(// , fih) + 2a 7 ^>(//, n h )] — 

x J \rriQ J 

x C D (m Q , » h ) f W du' S ?'^:!!?\ ; V = 2a r (/i, /x ) , (45) 

where Cj = Mh(1 — x)/x on the right-hand side. This relation allows us to derive the frag- 
mentation function Dq/ h from a primordial function Sq/h defined at a low renormalizat ion 
scale. 



4 Phenomenological results 

We now illustrate our results for the case of the 5-meson fragmentation function Di,/b(x, m&, fi) 
with the help of a phenomenological model. We use the function (35) with the parameter 
choices e B = 0.69 GeV, b = 2.9 (set 1) and e B = 0.82 GeV, b = 3.2 (set 2) to model the HQET 
fragmentation function Sq/ h (uj, fi ) at the scale /i = 1.5 GeV. We then use relation (45) to 
obtain the fragmentation function D),/b(x, to&, //) in the n;-flavor theory at a scale fi > /j,q. 
The evolution from the low scale to a higher scale is performed at next-to-next-to-leading 
order in renormalization-group improved perturbation theory, using the results collected in 
Appendix C. In order to test the stability of the perturbative approximation, we vary the 
hard matching scale Hh between m b /\/2 and \f2m\ ) . Our results are shown in Figure 3. We 
observe that the evolution to higher scales softens the fragmentation function and introduces 
a significant radiation tail. The dependence on the choice of fih is very weak. Depending on 
the set of input parameters the functional form of the fragmentation function can be varied. 
Parameter set 1 provides a good description of the B- meson shape function [22]. Parameter 
set 2 is chosen such that the resulting fragmentation function at the scale n = rrib resembles 
that obtained from a recent phenomenological fit to e + e~ data [52]. 

For practical purposes it does not matter whether a model for the fragmentation function is 
specified at a low or at a high scale. As long as one is interested in, say, £?-meson fragmentation 
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Figure 3: Scale dependence of the heavy-quark fragmentation function D b /B(x, m b , /•*) obtained 
from (45) using the model functions specified in the text. The overall normalization Nb is 
left free. In each figure, the three sets of curves refer to fj, — m& = 4.7 GeV (blue), \i = 3 GeV 
(red), and fj, — [Iq — 1.5 GeV (green). Within each set the curves show the residual dependence 
on the hard matching scale fih- The dotted line shows the phenomeno logical fit at fi = m b 
obtained in [52]. 



in e + e~ annihilation, one might as well stick with formula (15) for the cross section, evolve 
the partonic cross section down to a scale \i = m b (for details see [16]), and model the 
fragmentation function D b /B{x,fi = m b ) directly. At this stage, all perturbative logarithms 
of (1 — x) have already been extracted and absorbed into the partonic cross section. The 
factorization of the fragmentation function itself and its scale evolution become important 
only if one attempts to connect the fragmentation functions extracted in different processes 
with each other. In particular, treating the charm quark as a heavy quark, it is possible 
to relate the 6-quark fragmentation function of B mesons to the charm-quark fragmentation 
function of D mesons. In comparing the two functions we work for simplicity at next-to- 
leading order in a s , so that we do not need to worry about the decoupling of real and virtual 
charm-quark pairs. 

The resulting relation is simplest if we choose to evaluate both fragmentation functions at 
a common scale [i and ignore resummation effects. In this case, we obtain at leading power in 
AQCD/m C) 6 



D b/B (x,n) 



C D (m b ,n) M B 
C D (m c ,fi) M D 



D c m{y,n) 



M 



y 



D 



Md + Mb 1 -^ 



(46) 



The relation becomes more complicated when the two functions are evaluated at different 
scales and resummation effects are taken into account. We then obtain 



D b/B (x, (i = m b ) = U (m c , m b ) 



M B e~™ / Md\ v f 1 dy' D c/D {y\u. 



m. 



M D T{ri) \mj J y y' 2 (i=y _ W 

y y' 



l—T) 



(47) 



where renormalization-group effects are included in the quantities rj = 2aY(m b ,m c ) and 



U {m c , m b ) = —— r exp 

C D {m c ,m c ) 



2S(m c , m b ) + 2a 7 s (m c , m b ) 



(48) 
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Figure 4: Heavy-quark symmetry prediction for the charm-quark fragmentation function 
D c /d{^i^ — m c) (solid blue line) obtained from the inverse of relation (47) and the phe- 
nomenological fit result for the 6-quark fragmentation function Df,m(x,fi = m&) determined 
in [52] and overlaid as a dotted line. For comparison we show as red dashed lines the fit results 
for the charm-quark fragmentation functions of D + and D° mesons obtained in [54]. 

which we evaluate in 4-flavor QCD. Using = 4.7 GeV and m c = 1.5 GeV, the evolution 
functions evaluate to r\ w 0.322 and U(m c , nib) ~ 0.965 at next-to-leading order. 

At present the extraction of the charm-quark fragmentation functions into D mesons ap- 
pears to be affected by large uncertainties, at least as far as the overall normalization is 
concerned. For instance, the authors of [20] and [54] find a factor of almost 3 difference in the 
normalization of the fragmentation functions of D + and D° meson. A reason is, perhaps, that 
a significant fraction of D mesons is produced indirectly via decays of D* mesons, and these 
must be subtracted in order to obtain the sample of primary D mesons. From a theoretical 
perspective, a large difference between the fragmentation functions of D + and D° mesons 
would be very difficult to explain in QCD, as the two functions must coincide up to small 
electromagnetic and isospin-breaking corrections. We will therefore use the inverse of relation 
(47) to derive the charm-quark fragmentation function from the fragmentation function of 
b quarks. The corresponding relation is obtained by replacing m c <-> mj, and Mp. <-> Mb 
everywhere. Since the exponent rj = 2ar{m c ,mb) is negative in this case, it is necessary to 
regularize the integral using a star distribution (see eq. (47) in [15] for details). The evolution 
functions now evaluate to rj —0.322 and U(m c ,mb) ~ 0.715 at next-to-leading order. 

Relations (46) and (47) are model-independent consequences of heavy-quark symmetry [39], 
which are analogous to the relations between semileptonic form factors of D and B mesons 
obtained by Isgur and Wise [53]. Identical relations hold for other pairs of fragmentation 
functions, such as those for D* and B* or D s and B s mesons. In practice, we expect these 
relations to receive sizable power corrections, since the charm quark is not very heavy on the 
QCD scale. 

To illustrate relation (47) and its inverse we use the result for Db/si^, ^ = rrib) obtained in 
[52] from a fit to e + e~ data and derive from it a prediction for the charm-quark fragmentation 
function D c /d(x, [i = m c ). In Figure 4, we compare this prediction with the phenomenological 
fit results for the charm-quark fragmentation functions of D + and D° mesons obtained in 
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[54]. In the region of large x values our prediction comes close to the phenomenological fit 
for D° mesons. Note that for values of x below 0.79 the parameter Co exceeds 0.5 GeV for the 
case of charm, so that we would expect significant corrections from subleading fragmentation 
functions neglected in our analysis. 

5 Conclusions 

We have performed a systematic factorization analysis for the fragmentation functions of a 
hadron H containing a heavy quark Q. Short- and long-distance contributions associated with 
different momentum scales have been separated from each other using effective field-theory 
methods. In several aspects our approach goes beyond previous analyses of heavy-hadron 
fragmentation functions and introduces some novel features. In particular, we have shown that 
the popular decomposition of the heavy-quark fragmentation function into perturbative and 
nonperturbative components, in which the first component is identified with the perturbative 
fragmentation function Dq/q in (7), does not provide a proper factorization of short- and 
long-distance effects. 

We found it useful to integrate out heavy-quark pair production (virtual and real) by 
matching the set of fragmentation functions D a /n(x, toq, /i) defined in QCD with ni + 1 active 
flavors onto a single nonperturbative fragmentation function Dq/h{x, tuq, fi) defined in "par- 
tially quenched" QCD with ni active flavors, in which the heavy quark is still present but pair 
production of heavy quarks is forbidden. Besides the reduction to a single hadronic entity, 
which offers practical advantages for the modeling of fragmentation functions, in this step all 
calculable, short- distance dependence on x is extracted from the fragmentation function. We 
have stressed that most logarithms of (1 — x) in the expression for the perturbative fragmen- 
tation function Dq/q(x, mq, ji) are of long-distance origin and thus cannot be extracted and 
resummed in a meaningful way. We have completed this first matching step at two-loop order 
with the exception of a single coefficient function, Cq(x) in (12), for which the terms that are 
regular for x — > 1 have not yet been determined. To complete the calculation of this coefficient 
is an important task. 

In the second part of the paper we have considered the endpoint region, where (1 — x) 
is a small parameter, typically of order Aqcd/^q- In this case the residual dependence 
of the fragmentation function on the heavy-quark mass tuq, which remains after the first 
step, can be extracted by matching "partially quenched" QCD onto heavy-quark effective 
theory. We have derived evolution equations for the relevant quantities in the low-energy 
effective theory and solved them analytically. The resulting relation (45) allows us to control 
the scale dependence of the fragmentation function in the large-x region analytically, and 
to derive symmetry relations between the fragmentation functions for charm and bottom 
hadrons. Examples of such relations, which are model-independent consequences of QCD in 
the heavy-quark limit, have been discussed in Section 4. 

The results obtained in this paper provide the first step in a systematic analysis of frag- 
mentation of heavy hadrons at large x in e + e~ or pp collisions. In a forthcoming paper, we 
will perform the threshold resummation of large logarithms in the fragmentation cross section 
near x — > 1 directly in momentum space. This will set the basis for a precise determination 
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of the functions D b / H (x, m b , /i) and D c / H (x, m c , //) from fits to experimental data. Once these 
functions have been determined at a scale /i ~ tjiq, the entire set of fragmentation functions 
D a /n can be obtained using relation (10). These functions can then be evolved up to higher 
scales by solving the DGLAP evolution equations (5). 
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Appendix A: Two-loop matching coefficients 

We collect the two-loop expressions for the matching coefficients C a /Q defined in (10) and 
given in perturbative form in (11). The two-loop coefficients c a can be extracted from [8, 9]. 
Since the analytic expressions obtained in these two papers are lengthy, we will not reproduce 
them here but rather refer to the equations where they can be found. 

We begin with the two-loop coefficient c g describing the conversion of a massless gluon into 
a heavy quark Q. We obtain 



c g {x) = df{x) + -T F \x 2 + {I - x) 



hi 



2 V 



m, 



C F T F In \ \x 2 + (l 



X 



X 



1 — X 



In 



1 - x) 2 m 2 Q 



- 1 



(49) 



where d g 2 \x) is given in eq. (19) of [9]. The second term arises from the application of the 
matching relation (9) for the running coupling constant in the one-loop expression for the 
gluon fragmentation function, while the convolution term arises from the one-loop cross term 
in (10). 

The matching coefficients c a with a = q,q,Q for the conversion of a light quark, light anti- 
quark, or heavy anti-quark to a heavy quark start at two-loop order and are directly related 
to the corresponding perturbative fragmentation functions. We have 



c q {x) 



Cg(x) 



d?\x), 



x) = d ( ?\x) , 



where df\x) and (1q\x) are given in eqs. (54) and (56) of [8], respectively. 
We finally focus on the coefficient cq, for which we obtain 



(50) 



c Q (x) = C F T F {F ( 



(C F T F ) 
Q 



+ C F [Cf 



1 + x z 



In 



x 



(1 — x) 2 m 2 



(51) 
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(a) (b) (c) 



Figure 5: (a,b) Diagrams involving heavy-quark pair production, which yield the contribution 
AFq(x) in eq. (51). (c) A diagram with identical color structure, which does not involve pair 
production and hence does not contribute to Cq. 

where Fq f f (x) is given in eq. (62) of [8]. The term AFq(x) contains those parts of the 

coefficients Fq f and Fq a f in eq. (59) of [8] that are related to the Q — > QQQ fragmentation 
channel. They are given by the cuts of the first two diagrams shown in Figure 5. Due to the 
presence of other graphs with identical color structure, these contributions cannot be extracted 
without a dedicated calculation. 



Appendix B: Relations between distributions 

Using the usual definition of plus distributions together with the definition of the star distri- 
butions in (24), it is possible to derive general relations between the two types of distributions, 
which are valid for test functions depending on variables x or uj related by Cj = mg(l — x)/x. 
We obtain 



'!^M) M ^ 



1-x 

ln(l — x 
1 — x 



_ l n JL S(l - x) 

m Q 



dx 
x 



(52) 



+ m Q 



1 



1 — x 



iln 2 ^5(l-x)-i^-)^ 
2 vtiq 1 — x J x 



Expressions involving plus distributions of complicated functions of x can be simplified 
using the identity 



m 

l-x 



ln n fl - x) 



/(*) 



1 - x 



ln n (l - x) 



+ 



+ 



6(1 -x) f 
Jo 



W®-™ in-(l-x), 
1 — X 



(53) 

which is valid for integer n > 0. Applying this result, the one- loop expression for the per- 
turbative heavy-quark fragmentation function in the first equation in (7) can be rewritten 

as 



D Q / Q (x,mQ,fi) =6(l-x) 

C F a s 



1 + 



, In -^V + 2 

2tt V 2 rrin 



+ 



2ir 



ln 



1 — x \ m Q (l — x) 2 



(54) 
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Using our results (23) and (25), and expressing the star distributions in terms of plus distri- 
butions with the help of (52), we find 



g(x) C D (m Q , n) S q/q (uj, /i) x 



dx 



6(1 -x) 



2vr \2 m 2 n 



(55) 



CpOl s , . 
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- 1 



+ 



4 In a; 
1 — x 



Appendix C: Renormalizat ion-group functions 

The exact solutions (42) and (44) to the evolution equations can be evaluated by expanding 
the anomalous dimensions and the /9-function as perturbative series in the strong coupling. We 
work consistently at next-to-next-to-leading order (NNLO) in renormalization-group improved 
perturbation theory, keeping terms through order a 2 in the final expressions for the Sudakov 
exponent S and the functions and a 7 s. We define the expansion coefficients as 



r CUS p(a s ) = r ^ + r\ (^j 



4tt 



(3(a s ) = -2a s 



+ 



a, 



+ ... 



(56) 



and similarly for the other anomalous dimensions. In terms of these quantities, the NNLO 
expression for the renormalization-group function ar in (43) is given by 



«r(^,/i) = 



+ 




47T 

32vr 2 



+ . . . 



(57) 



The NNLO expression for the Sudakov exponent S in (43) has been given in eq. (94) of [15]. 
We do not reproduce this lengthy formula here. 

In the following we list the relevant expansion coefficients for the anomalous dimensions and 
the QCD /3-function, which are required at NNLO. All results refer to the MS renormalization 
scheme. For the convenience of the reader we also quote numerical values corresponding to 
rif = 4 massless quark flavors, as relevant for studies of 6-quark fragmentation. 

The expansion of the cusp anomalous dimension r cusp to two-loop order was obtained a 
long time ago [47] , while the three- loop coefficient has been calculated in [48] . The results are 



To 
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4C F 
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27 3 



+ C F T F n f (-f + 16C 3 )-^T|nJ 



429.507. 



(58) 



One also needs the four-loop cusp anomalous dimension, for which the Pade approximation 
T 3 ps 4313 (for rif — 4) was obtained in [55]. 

The anomalous dimension 7^ is know to three-loop order [48]. The expansion coefficients 

are 



lt = 3C F = 4, 

it = C 2 F (J - 2tt 2 + 24( 3 ) + C F C A ^ + 



17 22tt 2 



r t I 2 , ^ 
C^n, - + - 



43.8275 , 



7 | = C 3(| + 3 ^ + ^ + 68C3-^C3 -240 ( ; 



, ri ( 151 205- 2 247tt 4 
9 135~ 

20tt 2 116tt 4 272 



844 A 8vr 2 A " 
^-Cs + — C3 + I2OC5 



+ C'T^nj -46 + — — + 



9 



135 



+ C F C\ 



1657 2248tt 2 tt 4 1552 



36 



+ 



81 



18 



9 



C3 + 40C5 
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617.965. 



(59) 



The soft anomalous dimension 7 s is known to two-loop order [32, 42]. Its expansion 
coefficients are 



7c? = -2C F = -f , 



7? = C f 



110 tt 2 io \ [8 2tt 2 \ ' 

+ — - I8C3 ] C A + ( — + — ] T F n f 



27 18 



27 9 



-61.4198. 



(60) 



We use the simple estimate 7I ps |(7f) 2 /7g ~ —1768 for the three-loop coefficient, as this 
formula works remarkably well for the known three- loop coefficients of the quantities r cusp , 
7* and p. 

Finally, the expansion coefficients for the QCD /^-function to four-loop order are 
r 11 r 4 t 25 

y 3 7 = y ' 
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A = — C\ - — C A T F n f - ACpTpiif « 51.3333 , 

A = m cl + {2d - 2 f c F c A - If cj) r Fn/ + (f c F + f c„) r> ? 

w 406.352 , 

149753 /1078361 6508 , \ /50065 6472 \ 2 1093 , 

/3 3 = + 3564C 3 - + Cs )n f + + Cs n 2 f + n 3 f 

1 A 6 S V 162 27 S J } V 162 81 S J f 729 f 

^ 8035.19. (61) 
The value of @ 3 is taken from [56] and corresponds to A^ c = 3 and T F = \. 
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